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To understand orbital-angular-momentum contributions is becoming crucial for clarifying nucleon- 
spin issue in the parton level. Twist-two structure functions hi and hi for spin-one hadrons could 
probe orbital-angular-momentum effects, which reflect a different aspect from current studies for 
the spin-1/2 nucleon, since they should vanish if internal constituents are in the S state. These 
structure functions are related to tensor structure in spin-one hadrons. Studies of such tensor 
structure will open a new field of high-energy spin physics. The structure functions 6i and 62 are 
described by tensor-polarized quark and antiquark distributions StQ and Srq- Using HERMES data 
on the 61 structure function for the deuteron, we made an analysis of extracting the distributions 
Srq and Srq in a simple a;-dependent functional form. Optimum distributions are proposed for 
the tensor-polarized valence and antiquark distribution functions from the analysis. A finite tensor 
polarization is obtained for antiquarks if we impose a constraint that the first moments of tensor- 
polarized valence-quark distributions vanish. It is interesting to investigate a physics mechanism to 
create a finite tensor-polarized antiquark distribution. 

PACS numbers: 13.60.Hb, 13.88.-|-e 



Origin of nucleon spin has been investigated exten- 
sively after the EMC discovery that almost none of nu- 
cleon spin is carried by quarks lll| . Recent studies of po- 
larized parton distribution functions (PDFs) are found in 
Ref. \2\ . Although a gluon-spin contribution is not deter- 
mined accurately, orbital angular momenta are likely to 
be the crucial quantities in explaining the nucleon spin. 
Such contributions have been investigated recently by 
generalized parton distribution functions in lepton scat- 
tering Q and will be studied possibly at hadron facilities 
i- 

There are other quantities which are sensitive to the 
orbital angular momenta. For example, there are twist- 
two structure functions 61 and &2 in spin-one hadrons 
[3, Q . They could be related to the orbital angular mo- 
menta of internal constituents because they vanish if the 
constituents are in the S wave. Of course, they probe 
a different aspect of orbital-angular-momentum effects 
from the current ones for the nucleon because they are 
related to tensor-structure nature of spin-one hadrons. It 
is noteworthy that tensor structure is not understood at 
all in the parton level, which suggests that a new field 
of spin physics should be created by investigating the 
tensor-polarized structure functions. 

New polarized structure functions (61, 62, 63, and 64) 
were introduced in describing lepton deep inelastic scat- 
tering from a spin-one hadron [5|-l7| . A useful sum rule for 
the twist-two function 61 was proposed in Ref. [8], and 
it is partially used in this work. In conventional hadron 
models, such tensor structure arises due to the U-state 
admixture [a, [a, Q , pions [l3 , and shadowing effects [ll[ 
if the target is the deuteron. However, the tensor struc- 
ture would not be simply described by such conventional 



models at high energies according to our experience on 
the nucleon-spin issue. On the other hand, a theoret- 
ical formalism was developed in Ref. [12| to investi- 
gate the tensor-polarized distributions at hadron facilities 
by Drell-Yan processes with polarized deuteron. There 
are related theoretical studies such as new fragmenta- 
tion functions ^M, , generalized parton distributions ^ | , 
target mass corrections |lq. positivity constraints [l6|, 
lattice QCD estimate [I4I, projection operators of &i_4 
|18j for spin-one hadrons. The first measurement of the 
structure function 61 was done by the HERMES collab- 
oration in 2005 pLil]. 

The purpose of this research is to propose a simple 
parametrization for the tensor-polarized quark and anti- 
quark distribution functions by analyzing the HERMES 
data. It is intended to understand the current status of 
the tensor distributions. Obtained distributions could be 
used for comparing them with theoretical model calcula- 
tions and for proposing future experiments. 

The structure function bi is defined in the hadron ten- 
sor W^iz [Sll3- It is expressed in term of tensor-polarized 
distributions (^t? and Stq) as [y, §, U^ 



6i(x,g2) = Jj^e^ [STq^ix,Q^) + ^Tq^ix,Q^)] 



(1) 



where i indicates the flavor of a quark and e^ is the charge 
of the quark. The variables Q^ and x are defined by the 
momentum transfer q asQ^ = —q^ and x = Q^ / {2M mv) ^ 
where Mm and v are the nucleon mass and the energy 
transfer, respectively. Hereafter, the Q^ dependence is 
not explicitly written in the PDFs. In this work, the 
bi, Srqt, Srqi, and unpolarized PDFs are defined by the 
ones per nucleon for a nuclear target, namely they are 



divided by the factor of two if it is the deuteron. The 
tensor-polarized distribution StQ is defined by 



Srqiix) = g°(x) 



j+\x)+qr\x) 



(2) 



where q^ indicates an unpolarized-quark distribution in 
the hadron spin state A, and it is also defined the one per 
nucleon. Namely, Srq indicates an unpolarized-quark dis- 
tribution in a tensor-polarized spin-one hadron. It should 
be noted that the notation StQ is not the transversity dis- 
tribution, for which similar notations (6 or At) are used 
in nucleon-spin studies, throughout this article. A sum 
rule exists for bi in a parton model [8|: 



dx bi{x) 



— lim tFnU) 



0, 



(3) 



if the tensor-polarized antiquark distributions vanish 
J dx Sxqix) = 0. Here, FQ{t) is the electric quadrupole 
form factor of a spin-one hadron at the momentum 
squared t. 

We analyze the HERMES hi data for the deuteron. 
The tensor-polarized distributions are introduced as the 
unpolarized PDFs in the deuteron (D) multiplied by a 
common weight function dTw{x): 



Srq^ix) = 6tw{x) qf^{x), 
5Tqf{x) = Oq Stw{x) qf (x) 



(4) 



Namely, certain fractions of quark and antiquark distri- 
butions are tensor polarized and such probabilities are 
given by the function dTw{x) and an additional constant 
aq for antiquarks in comparison with the quark polariza- 
tion. The x dependence of StwIx) for antiquarks could 
be different in general from the one for quarks. However, 
it is not the stage of suggesting such a difference from ex- 
perimental measurements as it will become obvious later 
in this article. 

It is known that nuclear modifications are less than a 
few percent for the unpolarized PDFs in the deuteron 
[2l| . The tensor-polarized distributions cannot be de- 
termined within a few percent accuracy at this stage. 
Therefore, the nuclear modifications are neglected in q,p 
and qf . Then, the PDFs in the deuteron are written by 
a simple addition of proton and neutron contributions: 
qD = (gP + gj^)/2 and gf = (qf + q^)/2. Furthermore, 
isospin symmetry is assumed for relating the PDFs of 
the neutron to the ones of the proton: u„ = d, dn = m, 
Un = d, and d„ = u. Then, the tensor polarized distri- 
butions are 

6Tq^{x)=6TU^{x)=6Td^{x) 

Uy{x) + dy{x) 



= 5tw(x) 



STq'^ix) = (5tS^(x) = STd'^ix) = (5ts"(x) == Srs'^ix) 

, , , 2u{x) + 2d{x) + s{x) + s{x) .. 
= aqdTw(x) , (5) 



where flavor-symmetric tensor-polarized antiquark distri- 
butions are assumed. Tensor polarized heavy-quark dis- 
tributions are neglected in this work. To be precise, the 
distributions extend to x = 2 in the deuteron, whereas 
the kinematical limit is x = 1 for the nucleon. There- 
fore, the tensor-polarized distributions given in Eq. ([5]) 
cannot describe the region at 1 < x < 2. However, the 
PDFs are very small and it is not the stage to investigate 
the tensor distributions in such a large-x region. 

We analyze the data in the leading order (LO) of the 
running coupling constant a^. The structure function bi 
is then given by 



bi{x) = —STwix) [5{uy{x) + dy{x)} 

+4:aq{2u{x) + 2d{x) + s{x) + s{x)} ' 



(6) 



The unpolarized PDFs Uy(x), dv{x), • • •, s{x) could 
be taken from a recent global analysis, for example, by 
CTEQ J2^, GJR 0, or MSTW jM). In this work, the 
LO version of the MSTW parametrization is used. For 
the functional form of Stw{x), we note that there is a 
constraint from the sum rule in Eq. ([3]). The integrated 
tensor polarization, namely the first moment, should van- 
ish for the valence quarks. It indicates that there should 
be a node in the a;-dependent function, so that an appro- 
priate parametrization could be 



Stw{x) —ax (1 — xy{xo — x), 



(7) 



where a;o is the position where Srqvix) (and STqix)) van- 
ishes. If the first moments vanish for the valence-quark 
distributions, the constant xq is expressed by the other 
parameters as 



xo 



/q dxx''+^{l - xy{uy{x) + dy{x)} 
/g dxx''{l - xY{uv{x) + dy{x)} 



(8) 



The a, 6, c, and aq are the parameters to be determined 
from experimental measurements. 

The parametrization of Eq. ([7]) is motivated by the fol- 
lowing considerations. First, the parton model indicates 
the existence of a node as mentioned. Next, we expect to 
have smooth polynomial functional forms in the limits, 
a; — >■ and 1, as usual in unpolarized and longitudinally- 
polarized PDFs. In addition, the existence of the node 
and the functional form are, for example, supported by 
theoretical estimates of a convolution model, or so called 
binding model, where the Z?-state admixture gives rise to 
an x-distribution with a node in 6i including the tensor- 
polarized antiquark distributions [a, Ifij, ISj . 

From an analysis of the HERMES experimental data, 
the optimum function Stw{x) and ag are determined. It 
is obvious from the data that the 6i structure functions 
are not accurately measured to discuss scaling violation 
or even details of x dependence. Therefore, a simplifica- 
tion is made by ignoring the scaling violation. The Q^ 
value is fixed at Q^=2.5 GeV^, which is about the av- 
erage Q^ of the HERMES measurements, for calculating 



TABLE I: Determined parameters in Eqs. ([5| and ^. Q^ is taken (5^=2.5 GeV^ 



Analysis 



XVd.o.f. 



Xo 



Set 1 
Set 2 



2.83 
1.57 



0.378 
0.221 



± 0.212 
± 0.174 



0.0 (fixed) 
3.20 ± 2.75 



0.706 
0.648 



0.324 
0.342 



1.0 (fixed) 
1.0 (fixed) 



0.229 
0.221 



the unpolarized PDFs [23| in Eqs. ^ and 
two types of analyses: 



We made 



• Set 1 : Tensor-polarized antiquark distributions are 
terminated (ag=0). 

• Set 2: Finite tensor-polarized antiquark distribu- 
tions are allowed (a, is a parameter). 

Due to the lack of data at large x, the parameter c can- 
not be determined from the current data. We checked 
that the x^ value is not much affected by this parameter. 
Therefore, it is fixed at c = 1 in our analyses. 

The determined parameters are listed in Table ID It is 
obvious that the fit is not good enough (x^/d.o.f.=2.83) 
if the tensor-polarized antiquark distributions are termi- 
nated (set 1) by fixing the parameter as aq = 0. If we let 
this parameter free in the analysis, it is a reasonably suc- 
cessful one (x^/d.o.f.=1.57). It is interesting to find that 
the parameter aq is larger than one, which indicates that 
a significant tensor polarization exists in the antiquark 
distributions. 

Analysis results are shown in Fig. [1] in comparison 
with the HERMES experimental data. Only the data 
with Q^ > 1 GeV^ are included in the analyses. The 
set-1 curve is shown by the dashed curve, which does not 
agree with the data in the small- a; region (x < 0.1) with- 
out the antiquark polarization. The overall fit is success- 
ful only if the antiquark polarization is introduced (set 
2) as shown by the solid curve. Of course, the results 
depend on the assumed functional form including the as- 
sumption of using the common weight function Stw{x) 
for the quark and antiquark distributions. However, it 



would be reasonable as long as a smooth x-dependence 
is valid for the weight function ^^^(a;). 

The determined tensor-polarized distributions in Eq. 
([S]) are shown in Fig. [2]by using the parameters in Table 
m The amplitude is sightly larger for the valence-quark 
distribution of set 1 because the antiquark distributions 
are terminated by setting aq = 0. The antiquark dis- 
tribution StQ^ is shown by the dotted curve and it is 
mainly distributed in the region x < 0.1. It is interest- 
ing to find that a finite antiquark tensor polarization is 
needed to explain the HERMES data on 6i . If its effect 
on the bi sum rule is estimated, we obtain 



dx bi (x) 



— MmtFoit) 
24t^o ^^ ' 



1 

18 



dx [SSruix) 



+ 2STd{x) + Sts{x) + 6ts{x) ] 
= 0.0058. (9) 

The choice of parametrization of Eq. ([7]) for the an- 
tiquark distributions could affect the numerical result. 
However, as it is obvious from Figs. [T] and [U the anti- 
quarks contribute only at small a; (x < 0.1). As long as 
the function dTw{x) is a smooth function at x < 0.1, the 
result is not significantly changed. 

This work is the first attempt to parametrize the 
tensor polarized valence-quark and antiquark distribu- 
tions. Including the antiquark tensor polarization, we 
obtained much smaller x^/d.o.f. and it led to a finite 
sum as shown in Eq. ^. This is a new and inter- 
esting result which needs to be explained theoretically. 
The integral is compared with the HERMES results [l^ , 
J°-^^^dxbi{x) = [1.05 ± 0.34 (stat) ± 0.35 (sys)] x 10"^ 
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FIG. 1: (Color online) Comparison with HERMES data ^T^. 
The solid and dashed curves indicate theoretical results with 
(oLq 7^ 0) and without (a, = 0) tensor-polarized antiquark 
distributions. The open circle is the data at Q^ < 1 GeV'^. 
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FIG. 2: (Color online) Determined tensor-polarized distribu- 
tions. The dashed and solid curves are the valence-quark dis- 
tributions xStq^ in the deuteron for set 1 {og — 0) and set 2 
{aq ^ 0), respectively, and the dotted curve is the antiquark 
distribution xSrq^ of set 2. 



and £^^ dxbi{x) = [0.35±0.10 (stat)±0.18 (sys)] x lO^^ 
in the restricted range with Q^ > I GeV^. The integral 
of Eq. ([9]) is similar to the Gottfried sum [2^ 



X 



F^(x)] = i + - 



[u{x)-d{x)], (10) 



where the deviation from ^[u^{x) — dj, (a;)]/3 = 1/3 indi- 
cates flavor asymmetric antiquark distributions. In the 6i 
case, the finite sum J dxbi suggests that a finite tensor- 
polarized antiquark distribution should exist. 

It is obvious from Fig. [1] that much better measure- 
ments are needed to investigate the details of tensor- 
polarized distributions particularly at medium and large 
X (> 0.2). Such measurements could be possible, for 
example, at JLab (Thomas Jefferson National Accelera- 
tor Facility) by measuring &i and also at hadron facili- 
ties such as J-PARC (Japan Proton Accelerator Research 
Complex) Hi] GSI-FAIR (Gesellschaft fiir Schwerionen- 
forschung -Facility for Antiproton and Ion Research) |27| 
by Drell-Yan processes with polarized deuteron [l^l. In 
particular, the Drell-Yan processes are suitable for di- 



rectly finding the tensor-polarized antiquark distribu- 
tions in Eq. dH). 

Summary: In this work, optimum tensor-polarized 
quark and antiquark distributions are proposed from the 
analyses of HERMES data on hi for the deuteron. We 
found that a significant antiquark tensor polarization ex- 
ists if the overall tensor polarization vanishes for the 
valence quarks although such a result could depend on 
the assumed functional form. Further experimental mea- 
surements are needed for hi such as at JLab as well as 
Drell-Yan measurements with tensor-polarized deuteron 
at hadron facilities, J-PARC and GSI-FAIR. On the other 
hand, it is interesting to conjecture a possible physics 
mechanism to create a finite tensor-polarized antiquark 
distribution. 
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